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Abhijeet R. Sonawane-"^ and Prashant M. Gade^ 

Center for Modeling and Simulation, University of Pune, Pune-4-11 007, India 

PG Department of Physics, Rashtrasant Tukdoji Maharaj Nagpur University, Campus, Nagpur-4-40 033, 
India 

We investigate coupled circle maps in presence of feedback and explore various dynamical phases observed in 
this system of coupled high dimensional maps. We observe an interesting transition from localized chaos to 
spatiotemporal chaos. We study this transition as a dynamic phase transition. We observe that persistence 
acts as an excellent quantifier to describe this transition. Taking the location of the fixed point of circle map 
(which does not change with feedback) as a reference point, we compute number of sites which have been 
greater than (less than) the fixed point till time t. Though local dynamics is high-dimensional in this case, 
this definition of persistence which tracks a single variable is an excellent quantifier for this transition. In 
most cases, we also obtain a well defined persistence exponent at the critical point and observe conventional 
scaling as seen in second order phase transitions. This indicates that persistence could work as good order 
parameter for transitions from fully or partially arrested phase. We also give an explanation of gaps in 
eigenvalue spectrum of the Jacobian of localized state. 
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Coupled map lattices have been investigated ex- 
tensively in past few decades. The dynamics can 
show very novel features when local dynamics is 
higher dimensional. Using feedback to make the 
local dynamics higher dimensional makes it pos- 
sible to systematically increase dimensionality of 
local map without changing several other features 
such as location of fixed point. Investigation of 
phase diagram of such a system shows possibility 
of existence of partially arrested state in which 
certain sites are near the fixed point while others 
are not. We study the nonequilibrium phase tran- 
sition from such a state to a state of spatiotempo- 
ral chaos and propose that persistence works as 
an excellent order parameter for such transitions. 



I. INTRODUCTION 

Pattern formation in spatially extended systems has 
been an object of extensive study in past few decades. 
The reasons for the interest in pattern formation are not 
far to seek. Pattern formation happens in several natural 
systems and scientists are interested in understanding it. 
Examples could be flame propagation or occurrence of 
spirals in a chemical reactions or patterns on the skins of 
animals modeled by reaction-diffusion processes. Several 
systems like biological patterns^, charge density waves or 
Josephson Junction Arrays^"—, lasers^ have been studied 
extensively from the viewpoint of dynamics and pattern 
formation. The practical importance of understanding 
above systems can not be overemphasized. Partial differ- 
ential equations, coupled ordinary differential equations, 
coupled oscillator arrays and coupled maps have all been 
used to model different physical and biological systems 



and have uncovered interesting modifications of equilib- 
rium, bifurcations and stability properties of collective 
spatiotemporal states. These systems have been stud- 
ied extensively^ii. However, while we have arrived at a 
good understanding of low dimensional chaos in past two 
decades, and understood routes to chaos in several seem- 
ingly disparate systems in an unified manner, the same 
is not true for spatiotemporal systems. 

As one would expect, there is practical relevance to 
these models. Certain simplified models of spatiotem- 
poral dynamics have been studied extensively in recent 
past. We would like to make a special mention of cou- 
pled map lattices which have been studied extensively in 
the hope that the understanding we have developed in 
low-dimcnsional systems could be useful in understand- 
ing spatiotemporal dynamics. Majority of the studies 
are in coupled one-dimensional maps. At times, they 
have been successful in modeling certain patterns ob- 
served in spatiotemporal systems. For example, the 
phenomena of spatiotemporal intermittency in Rayleigh- 
Benard convection has been modeled by coupled map 
lattices'^'^. They have been used to model spiral waves 
in B-Z reaction^- and even in crystal growth^^. 

We feel that these studies could be helped considerably 
by quantitatively identifying different phases in the sys- 
tem and attempting to understand transitions between 
those. A lot of attention has been devoted to transi- 
tion to synchronization in these systems which is easy 
to identify and analyze. However, other spatiotemporal 
patterns are far less analyzed and understood. 

It is worth noting that in several spatiotemporal sys- 
tems the field is high dimensional. For example, for a 
chemical reaction, we would require information about 
concentrations of all reactants and products at all points 
to be able to simulate the system. There are several phe- 
nomena which are possible only when dynamics is higher 
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FIG. 1. We plot the phase diagram of coupled circle map with 
delay in two parameter space of coupling strength e and delay 
time r. The abbreviation SFP, FRP, FP-P/C and STC de- 
note Synchronous fixed point, Frozen random pattern. Fixed 
point-Periodic/Chaotic cycles and Spatiotemporal chaos re- 
spectively. In the inset, we magnify the phase boundary be- 
tween FP-P/C regime and STC regime. We can have only 
integer values of feedback time and the lines in phase dia- 
gram are guide to the eye. The black solid circles indicate 
points in phase space where spatiotemporal fixed point phase 



dimensional. On the other hand most of the existing 
studies in coupled map lattices are about coupled one- 
dimensional maps and very few involving high dimen- 
sional maps, for example coupled Chialvo's map^^, cou- 
pled Henon mapi^, Arnold's cat mapi^, modified Baker's 
map for Ising modeli^, etc. Here we try to study cou- 
pled higher dimensional maps. We try to systematically 
change the dimensionality of the map and study cou- 
pled maps with feedback where the feedback time could 
be varied. In coupled map lattice (CML)^- time evolu- 
tion is conflict between competing tendencies: the dif- 
fusive coupling between units tends to make the system 
homogeneous in space, while the chaotic map produces 
temporal inhomogeneity, due to sensitive dependence on 
initial condition. Such systems display a rich phase dia- 
gram and origin and stability of various phases in these 
systems are of theoretical and practical interest. The 
transitions between various phases in these systems could 
be viewed as dynamic phase transitions. Though these 
are clearly non-equilibrium systems, different techniques 
used in the study of equilibrium phase transitions have 
been applied to explore the variety of phenomenologically 
distinct behaviors. Such non-equilibrium systems range 
from growths of surfaces^i to traffic jams in vehicular 
traffic and granular flow a^^'^^ . For such analysis, there 
is a need to define an order parameter that will charac- 
terize the transition. However, not many dynamic phase 
transitions have been studied from this viewpoint^^— . 
Most of the studies are devoted to transition in models 
with so-called absorbing states. In the context of syn- 



chronization, there have been extensive studies to ascer- 
tain whether or not this transition is in the universality 
class of directed percolation^^. Transition to synchro- 
nization, locking in of the dynamics of several systems 
to a collective oscillation, is an important but minuscule 
part of overall spatially extended systems. Several non- 
synchronization transitions have been observed such as 
spiral state to spatiotemporal chaos^S, traveling wave to 
spatiotemporal chaos^. The spatiotemporal dynamics is 
far richer and other transitions deserve attention. 

As suggested by Pyragas, feedback in the dynamical 
systems can play a central role in controlling chaos by 
stabilizing UPO's embedded in chaotic attractor— . The 
scheme has been efficiently implemented in wide range 
of experimental applications. Feedback leading to syn- 
chronization is a well established fact supported by sev- 
eral researchers21. Though feedback has been an impor- 
tant aspect of studies in control and synchronization, not 
much has been said about its effect on different dynam- 
ical phases. It is then important to find the effects of 
feedback in terms of alterations in system dynamics that 
may lead to origin of different dynamic phases. 

In this work, we investigate coupled maps with feed- 
back. As mentioned before, coupled higher dimensional 
maps have not been studied adequately. Using feedback, 
we can systematically vary dimensionality of local map 
using feedback and study its effect on overall dynamics. 
We note that the feedback as defined in our model does 
not change the location of fixed point of system with- 
out feedback. We have used it extensively as a reference 
point. This allows us to study spatially extended higher 
dimensional systems with known properties. We investi- 
gate a model of coupled circle maps with feedback. The 
strength of feedback and the delay time are the param- 
eters in the system. We will present a general phase 
diagram of the system. In particular, we observe a novel 
transition from localized chaos to spatiotemporal chaos. 
By localized chaos, we mean that the laminar and turbu- 
lent phases coexist in the same space domain for weaker 
couplings. The reason to call it localized chaos is that 
one can see there are well defined boundaries of turbu- 
lent sites which do not mix with laminar ones. A similar 
pattern has been named spatial intermittency (SI) in pre- 
vious studies and some of its properties such as eigenvalue 
spectrum of Jacobian are found to be related to its dy- 
namical properties^SiSi. Interestingly, we find that per- 
sistence acts in effective way to characterize and quantify 
certain dynamical phases and transition between them in 
this system. Thus, persistence can be taken as excellent 
order parameter for certain phase transitions. 

The plan of paper is as follows. We introduce the 
model in the section II. In Section HI, we make a detailed 
analysis of one of the phases, namely localized chaos and 
calculate quantities such as turbulent length distribution. 
In section IV, we present our investigations in transition 
between localized chaos regime and spatiotemporal chaos 
and we present evidence that persistence acts as an ex- 
cellent order parameter to characterize the transition. In 
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section V, we demonstrate that some characteristic fea- 
tures in eigenvalue distributions of Jacobian can distin- 
guish between different dynamical phases. 

II. MODEL 

We define model of coupled circle maps with feedback 
as follows. We assign a continuous variable Xi{t) at each 
site i at time t where 1 < i < N. N is the number of 
lattice points. The evolution of Xi{t) is defined by, 

x,it + l) = il-e~5)fix^{t)) + 

+f{x,+im+df{x,{t-T)) (1) 

The parameter e is the coupling strength and the local 
dynamical update function f{x) is the sine-circle map, 

k 

f{x) ^ X + ui sin(27rx) (2) 

Where, k is the strength of nonlinearity and u is the 
winding number. The parameter S is strength of feedback 
and T is delay time. If we increase t, we can increase the 
dimension of map. In our simulations, we confine values 
of parameters to 6 ^ 0.1, fc = 1 and u = 0.068. The 
coupling topology considered here is diffusive i.e nearest - 
neighbor connections. For coupled circle map, we confine 
the dynamics in the interval [0, 1] using the following rule. 
If int[xi(i)] = m, Xi{t) = Xi{t) — to if Xi(t) > and 
Xi{t) = Xi{t) — m + 1 a Xi{t) < 0. We study the system 
with random initial conditions. The system is updated 
synchronously with periodic boundary condition. The 
feedback does not change the fixed point solution of local 
map f{x) without feedback. The fixed point solution for 
the local map f{x) is given by, x* = ^ sm^^{^jr-)- The 
parameter values considered here yield the value of fixed 
point as x* = 0.07026. This point acts as an excellent 
reference point in visualizing the dynamical behavior of 
system and also in symbolic dynamics of the system. We 
can define the site dynamics to be laminar if \xi {t)—x* \ < 
T] for small enough neighborhood rj — 0.001 and to be 
turbulent otherwise. 



III. SPATIOTEMPORAL DYNAMICS 

We examine the phenomenology of delay coupled circle 
map with feedback under variation of coupling strengths 
and delay time. Various spatiotemporal patterns are 
clearly evident from phase diagram Fig. [1] and the bi- 
furcation diagram shown in Fig. [51 Note that on y axis 
of Fig. [1] the time is discrete. The delay time t can only 
be some integer. The phase boundary of various regions 
can be seen. 

Increasing value of r is equivalent to increasing the di- 
mensionality of local dynamics. There are clearly four 
distinct behaviors shown by the system as seen from 
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FIG. 2. Bifurcation diagram of Spatial pattern at an instant 
of time (after transience), with respect to coupling strength e 
for delay time r — 1. 

phase diagram Fig. [TJ The first region in phase dia- 
gram labeled as SFP (Synchronous fixed point) where, 
for smaller coupling strengths e < 0.06 the entire sys- 
tem goes to local fixed point x*. This is an absorbing 
state and system settles in this state very quickly. We 
have checked that this state exists at least for t < 10. 
For higher coupling strengths. The system goes to spa- 
tiotemporal chaos as shown in phase diagram marked as 
STC. These different behaviors can also be visualized us- 
ing the bifurcation diagram Fig. [2] obtained for different 
coupling strengths for t — 1. This bifurcation diagram 
shows the state of all sites at a given instant of time. 
For very high coupling strengths, the system shows some 
stability islands for some values of r though they consti- 
tute a small part in phase diagram, we find them worth 
a mention. 

In the intermediate coupling strengths, the system 
shows a rather interesting behavior. We can see that 
some of the sites approach local fixed point and some of 
the sites take on different values which do not change in 
time. This behavior is shown in phase diagram marked as 
FRF which stand for frozen random pattern'^^. This re- 
gion is characterized by a spatial sequence of different do- 
mains whose size and periodicity may greatly vary. Upon 
further increase of coupling strengths, one enters in a re- 
gion where we can see coexistence of different dynamical 
regimes. Where some of the sites indeed go to fixed point, 
some of the sites show periodic/chaotic behavior. This 
is an interesting state where one can see that different 
dynamical phases do not mix with each other. The re- 
gions marked as FRF and FP-P /C in the phase diagram 
collectively represent the dynamics which we refer to as 
Localized chaos. Chaos is localized in certain dynamical 
regimes of turbulence and does not interfere with evolu- 
tion of laminar states. In localized chaos regime, chaotic 
domains and periodic domains are fixed in space. This is 
evident from spatial distribution pattern. The perturba- 
tions do not propagate in space. Essentially, the dynam- 
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ics is non-spreading and non-infective in nature, i.e. the 
spatially intermittent solutions have zero velocity com- 
ponents in spatial direction, and modes traveling along 
the lattice are absent. This can also be considered as a 
case of spatial intermittcney^*^ . The important result is 
that spatial intermittency is seen over a large parameter 
regime in presence of feedback. For each feedback time, 
we have a well defined region where we get the transition. 
These spatiotemporal patterns are clearly evident from 
space-time density plots below. Each figure represent a 
particular region in phase diagram. 

Fig. [3] shows space-time density plots of Xi^t for six 
choices of parameter values. The site index is plotted 
along the x-axis and the time evolves along the y-axis. 
As mentioned previously, laminar regions are blue (dark 
grey) in this representation. The colorbar shows differ- 
ent value of x{i) according to given scale. Fig. [Sja) 
gives the space-time density plot for parameter values 
where the system goes to spatiotemporal fixed point. Fig. 
[3Kb) shows the plot for non infective and non spread- 
ing behavior called as frozen random pattern. Fig. [SJc) 
and (d) shows coexistence of different dynamical regimes 
for which the dynamics essentially remains 'arrested'. 
FiglSfe) shows density plot for e = 0.24097 which is crit- 
ical point for the transition, we can see the spreading 
behavior just emerging in the system. Figl^Kf) depicts 
fully grown spatiotemporal chaos and here we can see 
that the state variable encompasses much larger part of 
available phase space. 

The number of turbulent sites increase as coupling 
strength increases. Note here that, all those sites which 
are not laminar are considered turbulent. We calculate 
distribution of turbulent phase length. It is observed that 
turbulent length distribution is governed by the exponen- 
tial law, 

T(0 = Aexp(-/30 (3) 

where, exponent /3 is the rate of decay and / is the the 
cluster size of turbulent sites. Fig. shows Cumula- 
tive distribution of turbulent lengths for coupling param- 
eter values e = 0.1,0.16,0.23,0.245. We can see that the 
distribution for smaller coupling strengths falls of quickly 
as opposed to higher value of coupling strength where 
the exponential decay is slower indicating large number 
of turbulent phases. Fig. S^b) shows how the exponent 
/? varies with respect to e. There is a sharp increase in 
value of beta near e ~ 0.24 . 



IV. PERSISTENCE AS AN ORDER PARAMETER 

As it can be seen the system displays a dynamical 
phase transition i.e. from localized chaos regime to spa- 
tiotemporal chaos. It is important to find a good quan- 
tifier to quantify this transition. This is a partially ar- 
rested state and in partially or fully arrested dynamical 
states, we believe that persistence acts as a good order 
parameter. There have been couple of previous studied 



in persistence in coupled map lattices by Menon et al. 
and by Gade et al.. Menon et al. found that it was 
a good indicator for a transition from fully synchronous 
state to spatiotemporal chaos while Gade et al. found 
that its dynamical behavior is able to characterize even 
traveling wave state which can be seen as a partially ar- 
rested state in moving frame of reference^^i^. It was 
observed by us in preliminary version of these studies^S 
that persistence acts as an good quantifier. Recent re- 
sults by Mahajan and Gade show that local persistence 
can also characterize the transition from clustered state 
to spatiotemporal chaos in small world networks^^. This 
definition of persistence concentrates on a single variable 
and one would wonders if it is adequate for the high di- 
mensional local dynamics we are studying in this case. 
However, we find that persistence acts as an excellent 
order parameter for distinguishing between partially and 
fully arrested states in this system. This demonstrates 
a generic utility of persistence as order parameter even 
in transitions from partially arrested phase. Thus, it is 
possible that persistence is a good characteristic for sys- 
tems in which dynamics is partially or fully localized. 
Conventionally, persistence in the context of stochastic 
processes has been defined as the probability P{t) that a 
stochastically fluctuating variable with a zero mean has 
not crossed a threshold value up to time t^- . For a spin 
system with discrete states, such as Ising or Potts spins, 
it is defined in terms of the probability that a given spin 
has not fiipped out of its initial state up to time t. The 
power- law tail of P{t), i.e. 

Pit) ^ l/^^ (4) 

defines the persistence exponent 9. The exponent 
9 is dimension-dependent, non-trivial exponent. For 
spatially-extended systems, the time evolution of the 
stochastic field at one lattice site is coupled to its neigh- 
bors. This makes the effective single-site evolution es- 
sentially non-markovian and renders the nontriviality to 
the exponent 9. This nature of local persistence prob- 
ability relates it to infinite point correlation function in 
time direction'^^ . 

A natural generalization of the idea to CML as sug- 
gested by Menon et al^, defines local persistence in 
terms of the probability that a local state variable a^j^ 
does not cross fixed point value x* up to time t. With 
this definition, we study persistence in the coupled map 
lattice defined through eq ([T]) , via simulations on systems 
with various lattice sizes. The persistence defined above 
is referred as local persistence in literature. In this work, 
we study only local persistence and do not study any 
other definition of persistence. In this paper, by persis- 
tence we mean local persistence only. We start with local 
Xi distributed randomly over interval [0, 1]. We find the 
normalized fraction of persistent sites, defined as sites for 
which Xi^t lias not crossed x* up to time t. Thus, sites for 
which the sign of {xi,t — x*) has not changed unto time 
t are persistent at time t. This persistence probability 
P{t) is averaged over an ensemble of over 2000 random 
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FIG. 3. (Color onlme)Space time density plot of coupled circle map lattice with feedback defined through eqs. (l)-(2), in a 
system of size L — 200, for parameters k — 1, S — 0.1, uj — 0.068 for r = 1 and (a)e = 0.05, (b) e — 0.1, (c)e — 0.16, (d) 
e = 0.23,(e) e = 0.24097, and (f) e = 0.245 



initial conditions. 

Fig. [5ja) shows plots of P(t) at three different values 
of coupling parameter e. It can be easily noted that for 
values below critical point P{t) saturates. For the value 
above critical point i.e. in chaotic domain, P{t) decays 
exponentially. Exactly at the critical point of transi- 
tion between localized chaos to spatiotemporal chaotic 
regimes, P{t) decays as a power law, which defines the 
persistence exponent 9 = 1.31. The power law decay 
can also be seen for other values of r-S . Fig. \5lh) 
shows persistence distribution for r = 4, 7 at the crit- 
ical points ec = 0.2386 and €c — 0.23847 respectively. 
The observed power-law behavior in 1-|-1 dimension sug- 
gests the further investigation of parameters involved in 
phenomenological scaling laws. The distance from criti- 
cality A = |e — ed and lateral system size L can be used 
as parameters to study scaling properties of the system. 
One would expect a asymptotic scaling law of the form 

p(t) ~^-^F(^i~^^A'^||^ (5) 

where, F is a scaling function and z — i'w/v^l = 1.62 
is the dynamical exponent. We consider two cases to 



demonstrate the validity of this scaling form numerically. 

As stated by Fuchs et. al.— , a finite 1-1-1 dimensional 
contact process at criticality reaches the absorbing state 
within finite time so that there is a finite probability for 
persistent sites to survive forever. Therefore, persistence 
probability saturates at a constant value. One can eas- 
ily carry this analogy to non spreading regime in coupled 
evolution of delayed circle maps in localized chaos regime. 
When the dynamics is such that no site is laminar, the 
persistence decays to zero. The scaling function at crit- 
icality A = for various system sizes L, i.e. plotting 
P{t)L^^ against t/L^, shows a excellent data collapse at 
longer times for z — 1.62 and 6 — 1.31 in Fig. |6l 

For a continuous phase transition below (above) the 
critical point, we expect the persistence probability to 
saturate (decay exponentially). According to above scal- 
ing form, the curves should collapse in both cases when 
P{t)A~^''\\ is plotted against tA'^ii. This is indeed the 
case as shown in Fig. [T] The simulations were carried 
out for lattice size L = 5000. The size is large enough 
so that finite size corrections could be neglected. Also 
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FIG. 4. (a) The distribution of turbulent phase lengths for 
coupling values e = 0.14,0.16,0.23,0.245 for r = 1. The 
fit to this distribution is exponential with exponent /3. The 
ordinate axis is presented in logarithmic scale, (b) This plot 
shows exponent /3 vs coupling strength e for t — 1. 



FIG. 5. (a) Persistence probability P{t) is plotted as a func- 
tion of time t on a logarithmic scale for fe = 1, w = 0.068 and 
T — 1. The lattice size is 5000. The curves shown are for 
e = 0.24040 (upper), e = 0.24097 (middle), and e = 0.24300 
(lower). The line with slop 1.31 is plotted for reference, (b) 
(Color online) Log P(t) vs Log t for r = 4, r = 7. 



ensemble averaging has been performed over 5000 differ- 
ent initial conditions to ensure convergence. Persistence 
acts as an excellent order parameter for transition from 
spreading to nonspreading regimes. This point can fur- 
ther be substantiated by similar power law decay seen for 
6 = case of spatial intermittency, as seen by Jabeen and 
Gupte^i. Figure m shows persistence probability behavior 
for parameter values of coupled circle map without feed- 
back where spatial intermittency is seen. The persistence 
exponent observed is 1.2. We would like to note that 
persistence exponent obtained in transition from laminar 



state to spatiotemporal chaos (shown to be in DP uni- 
versality class) is found to 1.5^ while in transition from 
laminar state to traveling wave state in negatively cou- 
pled circle maps, the exponent is 1—. Thus the exponent 
is clearly different from the exponents observed in pre- 
viously studied transitions. However as noted previously 
in literature, persistence is known to be least universal 
exponent among various dynamical exponents. 




FIG. 6. (a) The persistence probability P{t, L) is plotted 
against time t for different lattice sizes L. (b)The scaling 
function/(a;) = P{t,L)L^^ is plotted against dimensionless 
scaling variable x = tL~^ . The data for different L values 
were found to collapse for longer times for z = 1.62. 



FIG. 7. (a) The persistence probability P(t, L) is plotted 
against time t for different coupling strengths e. (b) The 
scaling function/(a;) = P(t)A~^''ll is plotted against dimen- 
sionless scaling variable x = f A"! I . A good data collapse is 
obtained for — 1.51. 



V. STABILITY ANALYSIS 



Computing Jacobian matrix of the system is useful 
from several vievifpoints. Lyapunov exponents are es- 
sentially related to eigenvalues of product of Jacobians 
over infinite period. For a synchronized or periodic solu- 
tions, one can get considerable simplifications by analyz- 
ing the Jacobian matrix for those states. The localized 
state has been previously observed in coupled circle maps 
by Jabeen and Gupte^i and has been named as spatial 
intermittency and they observed certain differences be- 
tween eigenvalue distributions of one-step Jacobian ma- 
trix. It has been observed that eigenvalue spectrum of 
the Jacobian of the system at any time T gives a good 
indicator of this phase. We show that the eigenvalues re- 
main real for our systems as well. This spectrum shows 
distinct features in the phases concerned. 




t 



FIG. 8. (a) Persistence probability P{t) is plotted as a func- 
tion of time t on a logarithmic scale for k — 1, uj — 0.062 
and e = 0.26634. The lattice size is 5000. The persistence 
exponent 9 = 1.2. 
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For T = 1, our equations can be put in the form. The Jacobian Jt at time t is given by, 

At Bn 
^ In On 

Xi{t + 1) = (1 - e - 5)f{x,{t)) + - [f{x,^i{t)) where Oat is an A^" x TV Null matrix with all entries 0. 

^ In is iV-dimensional identity matrix. The matrix Bn is 

+f{x,+i{t))] + Sf{y,{t)) defined as B{i,i) = Sf{yi{t)) and B{i,j) = 0, if i ^ j. 

yi(t + 1) = Xi(t) The matrix At is given by, 



At = 



/ej'ixiit)) e„/'(x2W) 

enf'{xi{t)) ej'{x2it)) e^fixsit)) 

enf{x2{t)) esf{x3{t)) 

\enf'ixi{t)) 



enf{xN{t))\ 


enf'{xi{t)) 
enfixN-iit)) eJ'ixNit)) J 



where, es = 1 — e — (5, e„ = | and f'{xi{t) = 1 — 
Kcos{2TTXi{t)). where, Xi{t) is the state variable. The 
diagonalization of Jt gives the 2N eigenvalues of the sta- 
bility matrix. The eigenvalues of the stability matrix 
were calculated for values of state variable Xi (t) for all i 
at one time step on discarding sufficiently long transients. 

Fig. inja) shows that the eigenvalue distribution for 
e = 0.2, where spatial intermittency or localized chaos 
is seen. This distribution has gap or zero probability 
regions in the eigenvalue spectrum as reported earlier by 
Jabeen et al?^ . On the other hand. Fig. IHJb) shows 
that there is no such gap for coupling values e = 0.3 
where fully grown spatiotemporal chaos is seen. Consider 
matrix J/ given by 



A't B'^ 
B'n On 



values. 



where B'{i,i) = ^5f'{y,{t)) and = 0, if i ^ j. 

The elements of the matrix A' are A'{i,i) = A{i,i) and 
A'{i,j) — \/ A{i, j)A{j, i) which can be obtained from the 
transformation A'{t) = C'AC'^^ where matrix C'{i,i) — 

One can verify that B''^ — B. Thus, it is easy to 
see that matrix J'{t) and matrix J{t) have same char- 
acteristic polynomial and hence they have same eigen- 
C" \ 

C B' I similarity transform which 

relates Jt and J/. Now the matrix J' {t) is symmetric and 
has real eigenvalues. Thus, the matrix J{t) also should 
have real eigenvalues. The assumption has been that 
f'{xi) > 0,f'{yi) > which is justified for fc < 1. The 
feedback strength 5 is already assumed to be positive in 
our work. Thus A{t) and A'{t) as well as Jt and J/ are 
related by similarity transformation and since A'{t) and 
J't are symmetric, A{t) and Jt should have real eigenval- 
ues. Thus we can look for gap in eigenvalue spectrum 
for these matrices and we observe the gap in our case 



for T = 1. Gupte and Jabeen already noticed the gap 
for case without feedback. We would also like to note 
that these are stretching transformations and eigenvec- 
tors which are localized in one frame of reference remain 
localized in other frame of reference as well. This obser- 
vation will be useful in later section. Thus it is possible to 
look for gap in the eigenvalue spectrum even in this case. 
Unfortunately the argument breaks down for r > 1—. 

Now let us try to understand the reason for gap in 
eigenvalue spectrum in A{t) in Jabeen and Gupte's case. 
We choose a diagonal part of matrix as an unperturbed 
matrix and all off-diagonal entries as a perturbation. Our 
choice for eigenvector for the diagonal part is a vector 
whose l*^ component is unity and all other components 
are zero. We define A[ — Af + A^, where diagonal en- 
tries of Af are given by Af{i,i) — At{i,i) and rest of 
its entries are zero. Let us consider an idealized case in 
which fc*'' and N*^ sites are turbulent and all other sites 
are laminar. Similar analysis can be carried out for more 
that two turbulent sites as long as they are not adjacent 
sites. 



a; 



... a 6fc 
. . . bk ak bk 
... bk a b 




\ 00 00 



bN \ 














a b^ 
bN aN I 



where a = (1 — e)f'{x*), b = e„/'(a;*) and at — {1 — 
e)f'{xi), bi — en\/ f'{xi)f'{x*), where i = k and N sites 
respectively. The matrix Af has eigenvalue a with (TV— 2) 
fold degeneracy and two nondegenerate eigenvalues given 



by ak and gn. The perturbation is At = A'^ 



Af. 



If all 
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FIG. 9. The eigenvalue distribution for (a) e = 0.2, (b) e = 0.3 is plotted for 50 different realizations for lattice size L = 1000 
witfi bin-width = 0.001. Gaps are seen in localized chaos eigenvalue distribution whereas the eigenvalue distribution for 
Spatiotemporal chaos region shows no such gaps. The inset in figure (a),(b) show the magnified region of respective plots 



sites are laminar except k^^ and iV*'* site (since we have 
periodic boundary conditions, by appropriate relabeling 
of sites, one can always bring the site in turbulent state 
to iV*'' site), the matrix Af ~ diag{a,a,ak ■ ■ ■ ,a,aN) , 
the perturbation will be. 



/060 0000 
6 6 









6fe 
6fe 6fe 
6fe 6 




V67vOOOOOO 



bN \ 









bN 
bN J 



We perform a simple degenerate state perturbation 
theory to our perturbation matrix A^, we can consider a 
matrix P of {N — 2) eigenvectors of degenerate eigenval- 
ues as, 



/lO 000 0\ 
1 



1 
00000000 
1 
1 



1 
VO 0/ 



it can be seen that 



fO 6 
6 
6 



P^APP = 



6 


















\ 











6 
6 






























6 

6 6 

6 6 

6 














6 

6 / 



This gives us for each {N — 2) fold degenerate eigenvec- 
tor, a (A^ — 2) X (A^ — 2) matrix of perturbation A^. We 
can see that the above matrix contains two blocks. The 
first block is of (fc — 1) x (fc — 1) size and the second block 
has size {N — fc — 1) x {N — k — 1). The eigenvalues of 
this matrix will be first order corrections. The perturbed 
eigenvalues after first order perturbation for degenerate 
states in this case will be. 



26cos(6'i) 



(6) 



wher( 



.36 



It: 



yi = 



- J (fc-i)-Hi 



for / = 1, 
for / = 1, 



,fc - 1 
,iN-k 



1) 



(7) 
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We have two non-degenerate eigenvectors, an unit vec- 
tor in TV*'* direction = [0, 0, . . . , 1]-^ and an unit 



vector in fc direction Vk — [0, 0, . . . 1, . . . , 0] 



Now 



'jJ^A^vn — v^A^Vk = 0. Thus the eigenvalues Uk and 
ajv are unperturbed after first order perturbation. Now 
if these values are far from a ± 26, we will observe a gap 
in the spectrum. 

Now let us generalize this argument to finite fraction 
/ of turbulent sites. The number of turbulent sites is 
m = fN and they are interspersed in the lattice. We also 
assume that no two turbulent sites are adjacent. In that 
case, we have = (1 — e)f'{xi), bt = enif'{xi)), where 
i = h,h,l3, ■■■,lm sites respectively, where Zi, Z2, ••, are 
indices of m turbulent sites and Im site is A^*'' site which 
can be obtained due to periodic boundary conditions. 
Again the matrix will have eigenvalue a with {N — m) 
fold degeneracy and m nondegenerate eigenvalues given 
by oi, a2, am- The matrix P becomes matrix of {N — 
m) eigenvectors of degenerate eigenvalues. The matrix of 
perturbation A^ changes accordingly to {N — m) x {N — 
m) . We can see that the matrix obtained from operation 
P^A^P contains m blocks. The square blocks will be of 
corresponding size, i. e. for li, the size will be (h — 1), 
for I2 it is {I2 — h — 1) and in the same way for m*'' 
block corresponding to turbulent site Z,„, the size will be 



{I. 



I, 



-1) 



1). The generalized case for m turbulent 



sites goes as. 







h 




llT 




h- 


1 






I3- 


2 






u- 


3 







for / = 1, .., h — 1 
for ^ = 1, .., ^2 — ^1 — 1 
for I = 1, ..,?3 - ^2 - 1 
for ^ = 1, .., — ^3 — 1 
for I ~ 1, .., — Im-i — 1 



(8) 



In any case the eigenvalues are between [a — 25, o -I- 26] . 
Now let us look at nondegenerate states. Arguing similar 
to the case of two site, after first order perturbation, the 
eigenvalues are a/j, aj^, . . . o;^. If these eigenvalues are 
not within the range a ± 26, we will observe a gap in the 
spectrum. Thus, we predict that eigenspectrum of Jt will 
be given by a band of eigenvalues between a-|-26 and a— 2b 
superposed with values en{f'{xj)), where j is a turbulent 
site. Localized eigenstates for A'{t) will also be localized 
for A{t) since they are related by stretching transforma- 
tion. If we make the histogram of eigenvalues as shown 
in figure [ini we find that this is indeed true. We can see 
histograms of eigenvalues computed numerically and the- 
oretically from the eqns. |n]and[5]obtained by pertubative 
analysis for 20 turbulent sites in the lattice of 100 ele- 
ments. The values for turbulent sites given by e„/'(a;j)) 
can also be seen away from the band [a — 26, a + 26]. For 
general case, there is a band of eigenvalues followed by m 
non-degenerate eigenstates which are far from the band. 
We checked that these eigenstates are localized and the 
corresponding eigenvectors are demonstrated in Fig. [11] 
for case of two turbulent sites. 

On the other hand, for spatiotemporal chaos, we have 
a symmetric random matrix A'^ which has nonzero en- 



TABLE I. Table showing Eigenvectors showing peaks at cor- 
responding two turbulent sites in otherwise steady state, here 
300 and 500. In second column larger eigenvalues are given 
which match with the values shown in last column. 



i 


Eigenvector 


Eigenvalue 




500 


1 


1.0861 


1.0817 


300 


2 


0.4177 


0.4132 




I numerical eigenvalues 
I theoretical eigenvalues 



a+2b 



0,2 0.4 0.6 0.8 



1 1.2 



FIG. 10. (Color online) Overlapping histograms of the cal- 
culated eigenvalues. In red (gray), histogram of eigenvalues 
obtained by theory of pertubative analysis. In blue(black), 
histogram of numerically obtained eigenvalues of matrix A(t) 
is given. The values for turbulent sites given by €nf'{xj)) 
can also be seen away from the band [a — 2b,a + 26] . The 
agreement between theory and numerics can be clearly seen 
for degenerate as well as non-degenerate eigenvalues. 



tries on diagonal, upper diagonal, lower diagonal and 
elements Af{l,N) and Af{N,l). In this case, ratio of 
width of band to system size tends to zero as iV — > 00. 
The eigenvalue spectrum of such a matrix tends to stan- 
dard Gaussian distribution—. It is worth noting here 
that, persistence acts as a good order parameter for this 
transition for all values of t including r = 0, whereas it 
is not clear that dynamic characterizer of gaps in eigen- 
values spectrum might work for higher value of t. 

Similar calculations can be carried out for our Jacobian 
Jj. The localized eigenvectors explain why localized pat- 
terns are stable. The eigenvectors having largest eigen- 
value are localized. Thus any perturbation from laminar 
state is likely to be growing in the direction of local- 
ized eigenvectors which does not disturb values at nearby 
sites. The above analysis of eigenvalue distribution for lo- 
calized chaos is restricted to demonstrate the analogy to 
spatial intermittency. The eigenvalue spectrum clearly 
differentiates infective to non-infective phase. The above 
analysis explains why the eigenvalue spectrum is different 
for spatiotemporal chaos and spatial intermittencyj^ 



11 











— IslEigenveclor 

— 2nd Eigemeclor 

















100 200 300 400 500 



FIG. 11. The localized eigenvectors showing peaks at 300 and 
500 for lattice size 500. 



VI. DISCUSSION 

In this paper, we have simulated coupled circle maps 
with feedback. This allows us to study coupled sys- 
tems where local dynamics is high dimensional. We give 
an overall phase diagram and discuss various dynamical 
phases. We also give an all-site bifurcation diagram and 
remark that there are phases which can not be found 
from bifurcation diagram and we need to study a de- 
tailed simulated space-time plot. In particular, we make 
a detailed study of a dynamical phase transition between 
the phase of localized chaos and spatiotemporally chaotic 
state. We find various characteristics such as distribution 
of domains of turbulent sites. Phase transition from dif- 
ferent dynamical phases is demonstrated. We make an 
attempt to analyze this transition in a manner in which 
equilibrium transitions are analyzed. In order that we are 
able to do so, it is necessary to define an order parameter. 
We find that local persistence acts as an excellent order 
parameter to characterize the phase transition. In addi- 
tion to this, conventional scaling laws pertaining to finite 
size simulations and offcritical scaling add credibility to 
use of persistence as quantifier of the transition. Inter- 
estingly, same definition of persistence used for without 
feedback case^ works for our case with feedback. We also 
calculate one step Jacobian of the system and find gaps in 
eigenvalue distribution of system which is dynamical sig- 
nature of non-infectious nature of dynamics in 'localized 
chaos' regime. The important effect of feedback is that 
the parameter window for which spatial intermittency is 
seen, gets enhanced. It should be noted that persistence 
acts as a good order parameter for this transition for all 
values of t including t — 0, whereas it is not clear that 
dynamic characterizer of gaps in eigenvalues spectrum 
might work for higher value of r. We get this transition 
for several values of r. We have shown that persistence 
acts as an excellent order parameter for this transition for 



all values of r. We also try to study gaps in the eigen- 
value spectrum which work well for t — and t = 1 . We 
go a step ahead in explaining the origin of such gaps in 
eigenvalues spectra of one step Jacobian by perturbation 
theory analysis. Analytical work also gives an idea of 
when above analysis may not work. 

There is a possibility to experimentally investigate dy- 
namical phases in spatiotemporal systems. Certain elec- 
trical oscillators are well simulated with circle map a^^i"^" 
and it is relatively easy to implement feedback mecha- 
nism in electric devices. If one can define a persistence 
in an analogous manner in differential equations, it is 
possible to find persistence in spatiotemporal systems ex- 
perimentally. 

In context of localized chaos, we would also like 
to point to an interesting work by Sethia et. al4i 
where they observe 'chimera' states, first observed by 
Kuramoto^. These states have spatially disconnected 
regions of synchronization separated by regions of de- 
coherence. Thus there is a strong analogy between 
'chimera' states and 'localized chaos'. Several works 
study 'chimera' under various modifications for identical 
oscillators 4^"— and non-identical oscillators^. Coupled 
evolution of systems with delay induces effectively non- 
local contributions in the evolution of individual system. 
Feedback propagates in space and time and affects the 
dynamics at some distant site. It could be of interest 
to find some general definition of persistence applicable 
to these systems and investigate transition to 'chimera' 
state as a dynamic phase transition. Further work in this 
direction is pursued. 
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